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grangians whose variational derivatives do not depend on second-order derivatives of the 
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Silva of the Regge-Teitelboim procedure that was used to define properly the mass in the 
classical formulation of Einstein's theory of gravity. Numerous applications reproduce 
standard results obtained by other secure but mostly specialized methods. As a novel 
application we calculate the Bondi energy loss in five dimensional gravity, based on the 
asymptotic solution given by Tanabe, Tanahashi and Shiromizu, and obtain, as expected, 
the same result. We also give the superpotential for Einstein-Gauss-Bonnet gravity and 
find the superpotential for Lovelock theories of gravity when the number of dimensions 
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1 Introduction 



(i) Hamiltonian formalism and the conservation of charges 

The way to evaluate total energy in the Hamiltonian formulation of general relativity [1] 
was put on a sound basis by Regge and Teitelboim |29| for localized source^] in spacetimes 
that are asymptotically flat. The correct Hamiltonian H has the following particularities: (1) 
The field equations are derived from a volume integral generator. (2) The generator is zero as 
a consequence of the equations, i.e. on-shell. (3) For any solution of Einstein's equations with 
given boundary conditions, H reduces to a surface integral on the boundary and with the 
boundary at spatial infinity the Hamiltonian is equal to the total mass-energy of the sources 
E including gravitational energjj^J 

In the covariant Hamiltonian formalism or covariant phase space as described by Julia 
and Silv^] [I{Jj the Hamiltonian appears in a slightly different light: (1) H is now associated 
with a timelike global vector field £, and is an integral over a spacelike hypersurface of a 
vector densitjj^] J. (2) The generator of the field equations is an integral over a spacelike 
hypersurface of a vector density current W, a linear combination of variational derivatives 
which is thus equal to zero on-shell. (3) The Hamiltonian on-shell reduces to a surface integral 
on the boundary S of the domain of integration V. The surface integral is generated by an 
anti-symmetric tensor density U defined on the boundary; this is the superpotential. The 
surface integral is equal to the energy E provided the timelike vector field becomes pure time 
translations at spatial flat infinity. 

We have thus the following sequence of equalities^] in covariant phase space: 

> = + dud^ => H= f PdVf, = [ W^dVf, + / U^dS^; 

Jv Jv Js 

on shell = , d M > = 0; and H = E = j> U^dS^. (1.1) 

In both cases E is the time component of a 4-vector in Minkowski spacetimJ^J If the 
timelike vector field £ is replaced by a spacelike vector field which turns (fast enough) into 
space translations along one of the spatial axes at infinity, the surface integral will be a 

1 Localized matter or fields that decrease faster than 1/r at flat spatial infinity. 

2 The method has been used in a variety of contexts, for a relatively recent application of the Regge- 
Teitelboim method see, for instance, Jamsin [14| . 

3 See also [18] which contains most references to original works on the subject. 
4 The density character is indicated by a hat " * " over a symbol. 

Greek indices go from to N — 1, N is the dimension of spacetime with signature — (N — 2). iV > 4. 
6 Various examples of applications can be found in [18] . 
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component of the linear momentum along that axis. The way to generate angular momentum 
components is now obvious. 

The formalism holds for fields other than gravity with other global or asymptotic symme- 
tries which are associated with other constants of motion or conserved quantities or what are 
generally called charge^} The whole machinery works also on a hypersurface at null infinity 
and allows one to calculate time-dependent quantities of physical relevance like the Bondi 
mass [5], the Sachs linear momentum |31j and the angular- momentum |21| with its usual 
ambiguities in four dimensions but apparently with no ambiguity in 5 dimensions |34j . 

(ii) The superpotential 

The superpotential U generates thus the conserved or non-conserved charges associated 
with the symmetries or Killing vectors of the metric on the boundary. 
There are two problems with the superpotential. 

(1) It must be calculated in Minkowski coordinates^] in spite of the fact that spherical 
coordinates are more convenient on a sphere, but in other coordinates the surface integral 
most often diverges. To avoid the latter one needs to introduce a background [20], [23] or 
counter-terms [H] , [27] or use a regularization procedure [21] or some other device [2] . 

(2) More importantly, the superpotential is not unique because the Lagrangian is not 
unique. 

We shall remove the coordinate difficulty by introducing a background metric g. The 

utility of backgrounds was suggested by Rosen [30] in the 1940's and often used since then 

in various publications, but we shall not take it as far as Rosen himself. The background 

metric g is defined as follows: take the metric g to be a solution of Einstein's equations that 

satisfies the boundary conditions and remove the source terms from the solution. Make them 

equal to zero. The result must be the metric far from the sources and thus on or near 

the boundary (at spatial infinity) g = g. Notice that the definition of the background is 

perfectly covariant. We shall never use the background as a "background" metric anywhere 

else but on the boundarjJ^J so we shall not have to define a mapping of our spacetime onto the 

7 The word charge without qualification may have different meaning for different authors. Here we used the 
word charge as it is understood, for instance, in Hollands, Ishibashi and Marolf |JTJ as constants of motions, or 
non constant if at null infinity, associated with asymptotic (and global) symmetries. Iyer and Wald [T3] define 
Noether charges. These are not quite the same and are more like some generalizations of Komar's covariant 
charges. 

8 See for instance Goldberg's |10| calculation of the Bondi Mass. 

9 In some cases, in particular in cosmological spacetimes, a background turns out to be useful [4j. 



3 



background. If spacetime is asymptotically flat, g is the Minkowski metric in any appropriate 
coordinates used for the solution of Einstein's equations. We then define the superpotential 
with respect to that background: U(g) — U(g). This is a covariant tensor which reduces 
to U(g) in Minkowski coordinates because in those coordinates U(g) = 0. The background 
solves the coordinate problem and allows us to calculate the total energy E in arbitrary 
coordinates at infinity, most usually in spherical coordinates. 

The superpotential non-uniqueness due to the non-uniqueness of the Lagrangian density is 



a different problem. It is best resolved by the method suggested by Regge and Teitelbom 10 
[29]. The variation of the Hamiltonian should have no boundary terms by analogy with 
classical mechanics. Thus the Hamiltonian H should contain a boundary term § UdS but its 
variation 5H should not. This will define <5U(g). Since the effective Hamiltonian generator 
is a linear combination of variational derivatives, c)U(g) thus defined is independent of any 
divergence added to the Lagrangian. Hence the superpotential is defined on the boundary up 
to an integration "constant". A complete unique superpotential is then obtained by taking 
the integration "constant" equal to — U(g). 

(Hi) Superpotentials, Affine Gravity and the Palatini formalism 

Silva [32] generalized the Regge- Teitelboim idea about the superpotential to the covariant 
Hamiltonian formulation of a large class of field theories, those in which the Lagrangians and 
their variational derivatives depend at most on first-order derivatives of the field components; 
moreover, the Lie derivatives of the fields with respect to a vector field £ depend at most on 
first-order derivatives of £ as well. Therefore, application to metric theories of gravity was 
done in affine-GL(iV;IR) gravity. The more familiar and practical Palatini formalism would 
not do because the Lie derivatives of the connection coefficients, the T's, contain second 
order-derivatives of 

The Silva covariant generalization of the Regge- Teitelboim method as well as subsequent 
developments |19] and applications were presented in terms of forms, within the "covariant 
symplectic" formulation. Our first intention here is to develop the formalism in component 
notations. 



Affine gravity introduces an extra set of N vectors Q a which constitute a local reference 



frame at each point of spacetime. This frame is arbitrary and the Lagrangian is invariant 



10 See also Henneaux and Teitelboim [T^] for anti-de Sitter backgrounds. 
11 Latin indices which enumerate the vectors vary also from to N — 1. 
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under any linear transformations 9 a (x) = A a b (x)O b (x). One motivation for this detour was 
to ensure that all Lie derivatives depend at most on first-order derivatives of Affine gravity 
is however not necessary to obtain the same results [22]. Therefore in a second step we shall 
review how to get the superpotential directly in the Palatini formalism and in component 
notations, in affine gravity and in the Palatini formulation. 

The formula for the superpotential we arrive at uses both affine gravity and Palatini 
formalism and as far as we are aware of is new. We also give various applications, some of 
which are new. 

(iv) Comment 

We want to emphasize that the superpotential exists only on the boundaries, which may be 
at spatial infinity or at null infinity. The closed boundary may be anywhere but in practice 
it should always be taken far away from the sources of gravity otherwise the value of the 
Hamiltonian is meaningless because gravitational energy is not localizable. 

2 Lagrangian variations in Palatini variables and in affine 
gravity 

In the Palatini formulation of gravity theories the fields are the metric of spacetime g with 
components g^ v and the symmetric connection^] T with components = T x Ufl ; these fields 
are independent. The curvature tensor does not have all the symmetries of the Riemann 
tensor, in particular d^T^ x ^ d u T pX . To avoid any confusion we shall denote the curvature 
tensor, like in Eisenhart [8], by B\ pa and keep the usual notation of the Riemann tensor 
Rr vpa for when the connection coefficients are Christoffel symbols. B v is antisymmetrical 
in the last two indices only, B\ = —B va . Additional symmetry properties, similar to 
those of the Riemann curvature tensor, are: 

B \u P *) = , 5 = and D * B V = D p B »° - D ° B »p where B ™ = B ^ * Bva ' 

(2.1) 

An occasional semi-column represents a covariant differentiation. The curvature tensor itself 

B\ pa = 2{d [p T\ ]u + T\ [p T^ ]v ). (2.2) 
12 Non symmetric connections are often considered also. 
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Now consider a Lagrangian scalar density C which depends on the metric components and 
the connection coefficients 

L = ^-gL^\B\ pa ). (2.3) 



vpa >• 

£ is a scalar and a hat over a symbol means it is multiplied by sj— g = — det(g^). The 
matter Lagrangian for isolated systems is not necessary in what follows. This implies that we 
take the boundary of spacetime outside of the matter, in practice at spatial or null infinity, 
or that fields tend to infinity faster than 1/r. 

Field equations are the variational derivatives equated to zero. The variation of £ for 
arbitrary variations of the field components 8g^ v and ST x pi/ may thus be written in the 
following form 

5C = *t^SsT + 4^ r V + ( 2 - 4 ) 

The content of £^ depends on the Lagrangian density which is defined up to a divergence 
and will not concern us here. The variational derivatives are 

*4, = 4 C - and C.^ = J^ ] (2.5) 

d stands for partial derivatives and 5 for variational derivatives. 

In affine metric-theories, the central element of the description is a set of N 2 vector 
field components 9 a ^ associated with a local basis in the local tangent spaces. There are 
N(N + l)/2 scalar products denoted by ^ ab : 

i ah = gTP^w (2-6) 

The symmetric connection T is replaced by a spin connection u which is a vector: 

= o\{d,e b x + T\ u e b ») = e\D,e b \ (2.7) 

Q b is the inverse matrix of a u . Reciprocally, 

T\ u = e\(6 a x u a vb -d v 6 b x ). (2.8) 

Since T is symmetrical in jjlv, the anti-symmetric part of the right hand side is of course equal 
to zero. In terms of and w, the curvature tensor becomes 

b\ pc = 2e a H h u {d [p u: a a]b + w ;/ ff]6 ) = e a H\B \ pa . (2.9) 

B a , depends only on u and its first-order derivatives. 
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The Lagrangian density £, a functional of g and T, see (2.2) and ( |2.3| ), may be expressed 
in terms 6 a , 7 a6 and u a ,. However, it will make things somewhat simplei) 13 l up to a point, in 



what follows if instead of 6 a , 7 afe and u a b we use a , g and u; a b . It is thus understood that 

t = L{g»\T\ v ) = £(<r,0%w%). (2.10) 

The variations of the Lagrangian will be written as usual in terms of variational derivatives 
which provide the field equations. Thus 



in which 



dC 



de a „ 



dg> 

The explicit form of T>^, like that of £ is not important. 



SC 



5co a 



(2.11) 



(2.12) 



lib 



We shall now go back from (|2.11|) to (|2.4|). To this effect we first replace duj 11 ^ in (|2.11) 



by its variations deduced from (2.7) in terms of SO and 5T: 



(2.13) 



The third term on the right hand side of (|2.11|) becomes, with minor obvious derivations by 
part in the square brackets 



ii 



t>* b 8u\ h = ( Cf b D 



v * 59 a u + 



/ fj. u b 



P x t>*0 h v 5!%. 



(2.14) 



This formula will take a nicer form if we introduce the following self explanatory new 
notations: 



a A £/ 6 = 4 ; similarly t 



fj,b. 



c 



fW 



9>4" = .1* and g^*L vX = *£^. (2.15) 



With these new notations we may rewrite (2.14) as follows 



(2.16) 



Inserting (2.16) back into (2.11), we find that 

5t = *t^Sgr + C X ^5Y\ V + + a x D v tr) ™\ + (2-17) 

13 Though not in the spirit of affine gravity theory! 

14 Partial derivatives equal covariant derivatives of scalars, vector densities and anti-symmetric tensor den- 
sities. 
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(2.17) looks like (2.4) except that in the Palatini formulation, there is no fl-dependenca^] 



Therefore, the factor of 56°"^ must be identically zero: 

• CJ 1 + 6 a x D u £f u = or • + D v t£ v = 0. (2.18) 
This identity will be valuable later on: it connects a variational derivative in affine gravity 



with the divergence of a variational derivative in the Palatini formalism. Equation (2.17) 
reduces thus to 

bt = *C^5g^ + C^5T\ V + 9^, (2.19) 
which is the same as hA. We notice that and = 6 b x eJ 1 t b av depend only on g 



and B. Affine gravity and the Palatini formalism of Lagrangians that depend only on the 
metric and on the curvature tensor but not on their derivatives lead to equivalent equations 
as expected. The variational derivatives are identical when the vector fields 9 are in the 



directions of the local axis and 6 a „ = 5°;. 

A 1 A 1 



3 Symmetries and conserved currents in affine gravity 



(i) Diffeomorphism invariance in affine gravity 



Conservation laws related to the arbitrariness of the coordinates are normally derived 
following Noether's procedure. One considers the Lie derivative of the Lagrangian density 
5lC with respect to an arbitrary vector field £. The Lie derivative of a scalar density is 
a divergence 5l£ = d^{C^). In (2.11) we replace the variations of the field components 
5g^ v ,59 a ^ and boj a ^ b by their Lie derivatives as well. These do not contain higher than first- 
order derivatives of £: 



(D x g^ X -2g x ^D x ^ , 5 L 9\ = (D x 6\) £ A + 6 a M 



5 L u\ h = {D x u\ b )e + u\ b D^ 



(3-1) 



We may thus rewrite (2.11) with Lie derivatives as follows: 



(3.2) 



Suppose we take variations of the fields that are all zero on the boundary and integrate over spacetime. 
The volume integral should not contain any SO terms since 80 is arbitrary. 



s 



From (3.1) it is easily seen that the divergence is of a vector density linear in £ and D£, say 



of the form 



- V u = *J^ X + *U^D^ X . 



(3.3) 



The left-hand side is equally linear in £ and D£. Thus (3.2) may be rewritten as follows: 



(3.4) 



in which 



A x = *C^D x g^ + . CfD^ + lfD^ , A" A = -2*1* + .t* + Lfco \ b . (3.5) 

It is interesting to notice and easy to verify that — 2*£^ + *£>J* is, in affine gravity, the 
variational derivative with respect to 9 a ^ times 9 a x : 



r a* 

A 



-2*£^ + = ^-0° A 

A A QQa A 



(3.6) 



~f,U> 



Now in (3.4) replace^] £(x) by £(x)e(x); this can always be done because £ is arbitrary. Then 



expand in terms of derivatives of e(x). The left-hand side of (3.4) looks then like this: 



(A A £ A + Ay^ A ) e + (A^ A )ZV. 



(3.7) 



On the other hand, the right hand side is now of the following form: 

O^re + *U v »d„e) = (8 u *r)e + (*> + d^U^d^e + {tU^D^e. (3. 



Thus, since (3.7) and (3.8) are equal 



(A A £ A + A\D^ x )e + (A^ A )ZV = (d v *J")e + (*> + d v *U^)d^ + (*U^)D^e. (3.9) 
e is arbitrary and we may identify the factors of e, d^e and of D^ v e of both sides. There 



follows a set of "cascade identities" [16]. The factors of e reproduce of course identity (3.4). 
The factor of D^e does not appear on the left-hand side ; thus *U U>M must be antisymmetric: 



(3.10) 



Taking this into account, the identification of the factors of d^e leads to the following identity: 
* = *W» + faU'"' in which *W" = A M A e A = + Cfuo a Xb )i x . (3.11) 



3 Like in Julia and Silva |16j . I like this clever trick (JK)! 



The vector is conserved, i.e. divergenceless, on- shelf] (4" = £/ 6 = 0). The value of 
the integral of on-shell over a hypersurface V is equal to the integral on the boundary 
S of V of *U >1U . However the "generator" is not manifestly covariant; it depends on 

the choice of u. In particular in the Palatini formalism 6 a ^ = 5® and 6^9^uj a ^ b = P^,. Thus 
* JP, which has much of the properties of a generator of conserved charges, is unfortunately 
not covariant. The way around this problem which we now describe is given in |19j . 

(ii) Freedom of choice of the local frame 

In affine gravity theories the 6 a, s are N arbitrary vectors which play the role of local 
tangent frames but the theories are invariant for arbitrary linear changes of basis 9 a (x) —> 
A a b (x)O b (x). Therefore we can play the same game with this invariance as we played with 
the diffeomorphism invariance. 

Let X a b characterize an infinitesimal transformation denoted by 5\. Thus, 

S A e a ll = X a b b lt implying that <5 A u;% = -d^X\ - u\ c X c b + u;%A a c . (3.12) 
The Lagrangian density and the metric components are invariant under such transformations: 

5 A C = , W = 0. (3.13) 



Thus if 5 is replaced by 5\ in (2.11), we obtain, after some rearrangement of factors, 



(4"^ - 4**% + cr^\ c )x a b - 4%n = d u (**r a b \% +**u^%\%), (3.i4) 



The V 1 in (2.11) is linear in X a b and its derivative and therefore the right-hand side in (3.14) 
is of the form written above. 

We now play the same game as before by replacing X a b (x) with X a b (x)e(x) and expanding 
in terms of e, d^e and D^e, identifying the corresponding factors on both sides. The factors 
of e reproduce of course (3.14). The factor of D^e tells us that 



**U^ b X a b = **U Vfl = -**U^ U . (3.15) 



The equality of the factors of d^e leads to an equality similar to (3.11): 



* 



*J^ b X a b = **J^ = -CfX a b + d v {**U^). (3.16) 



7 With sources present the equations are of course different, except on the boundary. 
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Here is another non-covariant divergenceless vector on-shell with the properties of 



see (1.1). It is divergenceless but not covariant. If we now sum *J of (3.11) and **J just 



written, we obtain a new conserved current on shell with the same properties, which can be 
written as follows: 



> = + = c^^ + tf (u a vb i v -\ a b ) + d u U^ with = *U^ + **U^. (3.17) 



(Hi) A covariant conserved current 

The final step in [18] entails the breaking of the local frame arbitrariness. The 9 a, s are 



kept fixed like they would in a Palatini formulation by taking 



5 L e^ + 5 A e a fl = o , d,e\ = o and e\ = 8i => \ a b = -e\e b ^e , e a x e b v ^ b = r\ v . 

(3.18) 



Substituting the last two equalities in the second term of (3.17) on the right-hand side gives: 

4/^ v A + trd»e = Ar r u A + = c^dj*. (3.19) 

So, J is now of the form 

j» = W» + dyU^ with W» = £/£ A + t£ v D v t x . (3.20) 

The covariant is conserved on-shell (W^ = 0) and the charge 

Q = j> U^dS^. (3.21) 

Q is, of course, not yet well defined because the superpotential is not well-defined. We come 
back to this problem with an answer below. But we shall first show that the same current is 
obtained by applying Noether's prescription to the same Lagrangian density in the Palatini 
formalism with g and T as independent fields without introducing an additional set of N 
vectors 0. We shall follow 122! in this matter. 



3 Or, if one prefers, in an Einstein-Cartan formalism — ej) 
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4 Symmetries and conserved currents in Palatini formalism 



Conservation laws related to diffeomorphism invariance are derived in the same way as in 
affine gravity. The difference lies in an extra derivative of The Lie derivative 5lC = d^C^) 



as before. The Lie derivative of g^ v is the same as in (3.1 ). However the Lie derivative of the 
r's contains second order derivatives and is a tensor: 



S L T' 



[IV 



- R (jiv) V £ V + D (jw)€ • 



(4.1) 



The result is that instead of an identity like (3.4) there is one additional term on both sides 



in Davi X - Thus after replacing 6 by 5l in (2.4), using (3.1) and (4.1), we obtain an identity 



of the following form: 



k x e + k d ^ x + t^D^e = 9u (r x e + •u x ni D,e + v x vpa D {pa) e) (4.2 



in which 



k x = *t tlv D xS r-t^W ' , A/ = -2*£ 



-x ■ 



(4.3) 



Now in (4.2) we replace £(x) by £(x)e(x) and, as before, we expand in terms of e, c^e, D^ u e 
and D XfjLU e. On the left-hand side we have 



(A A £ A + A/£> M e A + Cjr D n»t)e + ( A xt + 2Cj^D v ^)D^ + (C^^D^e. (4.4) 
On the right hand side we have however this: 

{8j u )e + (> + d v mU^d^e + {•U v » + D x V x ^)D, w e + V XfM/ D Xfn/ e. (4.5) 



The factors of e and its derivatives on the left-hand side (4.4) are identical to the factors 



of those derivatives on the right-hand side (4.5). The identity "left-hand side (4.4)= right 



hand side (4.5)" holds in arbitrary coordinates. We may take coordinates in which at any 



particular point the metric is Minkowski and the Vs are equal to zero. At that point the 
last term in (4.5) becomes V Xflu d x ^ u e. It must be zero because there is no such term on the 



left-hand side (4.4). This implies that V x ^ u , totally symmetrized in all its indices, equals 
zero or that 



y{\»v) = Uy^v + + V uX ^) = and, see (Ob, V x ^ u = V Xv » 



(4.6) 
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Any V that has these properties has also the following property which is easily proven: 



With this equality (4.7), we may rewrite the right-hand side in (4.5) as follows 



(4.7) 



(4.8) 



The left-hand side (4.4) is thus also identical to the right-hand side in the form (4.8) for any 



values of e: the factors of e and its derivatives in (4.4) are identical to the corresponding 



factors in (4.8) 



The factors of e are obtained by setting e = 1 which reproduces thus exactly (4.2) from 



which we learn nothing new. The two remaining identities obtained by identifying the factors 
of <9 M e and D^e, are respectively 



-2 *C^ X + 2C^D U £ X = > + D U »U U ^ - lV^ v R^ 
£ = m uW + D X V X ^. 



VpCJl 



(4.9) 
(4.10) 



We shall decompose •U Ufl appearing in (4.9) into a symmetric •U^ u ^ and an anti-symmetric 



part •U^ Uf ^ in [iv. The symmetric part will be replaced by its value given in (4.10), i.e. 

• (j(^) = t£ v i x - D X V X » U . 



The right-hand side of (4.9) then becomes 



vpu ■ 



However, the last two terms combine to give 



D x V Xflu - IV^R^pa = \d v (D X V^ X ), 



so that the right-hand side of (4.9), i.e. (4.12), is now of the form 



> - djj» v + £ A ^£ A in which = -U vtx = •U^ v] - \D X V^ X . 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



Since the left-hand side of (4.9) is equal to (4.12), we have now 



> = W£ + djjW in which W» = (-2 * - Dxtj*)? + t£ v D v ? 



(4.15) 
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This expression is obviously the same as (3.20) in affine gravity since, if we take into account 

(4.16) 



the identity (2.18), 



which is indeed (3.20). 



The passage from affine gravity to the Palatini formalism is implicit in the work of Julia 
and Silva [18] . One difference is that the whole layout is here explicit and in tensorial 
notations. 



5 Superpotentials and charges in Lovelock gravity theories 



(i) Hamiltonian covariant formalism 



With an expression for the current conserved on-shell (4.15) we may write the covariant 



Hamiltonian H. We keep the letter H for in the end we are mostly interested in the total 
mass energy. However, at this stage we shall not specialize £ to an asymptotic timelike Killing 
vector. Thus, in the Palatini formalism, 



H 



U^dS, 



U^dS, 



flU- 



(5.1) 



'V JS JV x ' JS 

We can equally use the field components of affine gravity, before breaking the gauge invari- 
ance. The calculation is easily seen to lead to the same result once gauge invariance is broken. 



From (|3.6|) we know that is an ordinary derivative which depends on g and B only. On 



the other hand, in Lovelock gravity, C{~ are variational derivatives, see (2.5), which contain 



not only g and B but are linear and homogeneous in derivatives of the metric components 
D\g^ v but contain no derivatives of B. The covariant Hamiltonian equations are obtained 
from the variation of H which reduces thus to: 



5H 



fig 



pa 



5T X 



pa 



pa 



dV^ + 



{d v9 P°) 



BP ^P 1 ! 



d (8 V T- 



pa 



pa) 



dS^ 
(5.2) 



It is not at all obvious that the term in front of 5U^ U in the brackets on the right hand side is 
antisymmetrical in [lis, however it is. The prove goes as follows: consider the left hand side 



of (3.4) which is a divergence. Its variational derivative with respect to any field component, 



say y A , is thus equal to zero. This holds for any £. Replace £ by e(x)£. It holds also for any 
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e(x). The factor of D pu e must thus be equal to zero. This factor is dW^ / d{d v \y A \. Thus 
the factors in the bracket on the right hand side are anti-symmetric in \xv. 

(ii) Dirichlet boundary conditions and the superpotential 

Following Julia and Silva the functional equation for the superpotential is given by equat- 
ing the second bracket to zero: 

+ ^_„. (5.3) 



pa 

pa J 



(d v g pa ) d{d v T\ 
Superpotentials must satisfy all boundary conditions. The most common conditions in Ein- 
stein's theory of gravity are the Dirichlet boundary conditions g pa \s = 9 pa , or equivalently 
$9 pa \s = 0- We shall also identify the T's on the boundary with Christoffel symbols which 
is equivalent to taking D\g pa \s = or 8(D\g p<T )\s = 0. As a consequence on the boundary 
B\pv\s = R^apvi the Riemann tensor of the background, and 5B\^ v \s = 0. Finally the 
field equations themselves must remain satisfied. Thus, for isolated sources the boundary 
conditions are: 

5gP°\s = 5(D x g p °)\s = SB^s = 5C P = 5C^ = 0. (5.4) 



Equations (|5.3|) reduce thus to 

-5T x pa + 5U pu = 0. (5.5) 



The factors of ST X in (5.5) are entirely defined on the boundary. As far as the differential 



d (d v T x 



pa 

pa J 



equation for U is concerned, we may thus integrate equation (5.5), and obtain an expression 



linear in T x pa . The "constant of integration" will be minus the same expression with p(T 
replaced by T X pa for reasons explained in the Introduction. The formula for the superpoten- 
tial is thus 

^ = -^P) A V A^ = rV-rV (5.6) 



Finally, since £ is a function of g and B, we may also rewrite (5.5) like this: 



^ = 2^§£^AV (5.7) 

u pav 

This formula for the covariant superpotential is a blend of two derivatives: one with respect 
to at 7 and u> constant in affine gravity giving the generator C p tp and one with respect 
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to B x pau at g constant in the Palatini representation. As far as we know this result is new, 
simple and compact. The corresponding charge^ Q is 



Q = f U^dS,* (5.8) 



in which U^ u is defined in (5.7). 

6 Some examples 

(i) Charges in N dimensional Einstein's theory of gravitation 

This simple example will show how the machinery works. The Lagrangian is 

A = — B = —g^B, lv = —Mar ~ 5Vg pu )B x = —8!g ffV B x u r ntT ,. (6.1) 
2k 2k m 4k Ay xy ' vpa 2k xy v w v 1 

In four dimensions we take as usual for k Einstein's coupling constant k = In N 



dimensions 

2Sn-2Gn 



(6.2) 



where Sjy-2 is the surface of a unit sphere of dimension (N — 2) and Gat is a gravity coupling 
constant. In affine gravity, the field components are "y ab , 6® and oj a pb in terms of which, see 



(2.6) and (2.9): 



A = ^IW (V^ft) where 9 = det W ' T = det ^) and B< J" = 

(6.3) 

The term in parenthesis is independent of the #'s. The calculation of the derivative with 
respect to a ^ is thus very simple. One has, however, to be careful about symmetrizations 
and anti-symmetrizations of indices. The quantity that interests us, the generator appearing 



in (5.6), is here xC^^P which in a condensed form may be written as 

i4/r = v- (es l x9 v]p + r<ft M]p + e^g u]p ) # V ( 6 - 4 ) 



2 k 

From this follows that the superpotential Ui is 

UF = ^ (&9 U]P ^ P , + &g ap ^l + i u 9 p[p ^%) ■ (6.5) 



19 The superpotential given in j22] is not the same as the one given in (5.7 1. The reason is that we used 
an equality (3.20) (in [55]) which does not correspond to a covariant generalization of the Regge-Teitelboim 
procedure. It is rather equivalent to what amounts in Julia and Silva 48] to take X 7^ in their equation (10). 
The superpotential we found in [22] gives however correctly the charges for backgrounds that are maximally 
symmetrical which is the most common case and which explain that all the applications described in [22| are 
correct. The claim of the paper that it generalizes the Regge Teitelboim method along the line described by 
Julia and Silva is, however, wrong. 
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Notice the circular symmetry of indices fjL,a,v in (6.4) and [jl,v,o~ in (6.5). In exactly this 
form, but in Minkowski coordinates on a flat background and for £^ = {1, 0, 0, 0}, the super- 
potential for the total mass energy in Einstein's theory in 4 dimensions was published by von 
Freud [35] in 1939. 



A more standard form is obtained by noticing that 



i u g p ^A v \ a = - £>&*|"] (6.6) 

and then writing like this: 

= -L (dH^ - dH v] ) + ir&k v] with k v = g vp A% a - g pa A v pa . (6.7) 



In this form we recognize the KBL [23] superpotentia 20 The first term j^D^^ is half 
the Komar |25| superpotential. For applications formula (6.5) is simpler. Various physical 
properties including standard results for charges derived from the KBL superpotential have 
been reviewed in |18| . It gives among other things the Bondi mass [5T] for radiating systems 
in four dimensions. The KBL superpotential is thus valid at both spatial and null infinity 
and on any background that satisfies the boundary conditions. Here is a new application of 
the superpotential at null infinity in five dimensions. We have calculated the mass loss using 
the asymptotic structure of the solution of Einstein's equations given by Tanabe, Tanahashi 
and Shiromizu |33j . see also |34| . in "Bondi coordinates" and obtained their formula for the 
mass los9, 21 [ Some details of the calculations are given in Appendix A. 



(ii) Superpotential of Gauss-Bonnet Gravity in N > 4 dimensions 

Einstein-Gauss-Bonnet theories of gravity are usually considered in combination with 
Einstein's Lagrangian [see sub-section (iii) below]. Since we obtained already Ui we shall 
consider U2 derived from the Lagrangian density £2 which we shall write in a form similar 
to (6.3) the most convenient to calculate the generator iCrf'tp and which we took from 



Lovelock's paper of 1971 



The terms within parentheses depends on 7's and lJs and not on 0. The 0's are all in front 



of the parenthesis. Equation (6.8) provides the generator v^r^t? through ordinary partial 



20 On a flat background this expression was given in |20| . 

21 Adjusted to our coupling constant in front of the Lagrangian and our units. 
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derivatives of (6.8) with respect to the 0's. The result is worth rewriting in terms of g's 



and B's, since the superpotential U2 is also obtained by ordinary partial derivatives of the 
generator 2A/ t £' ? with respect B's at fixed g's: 



(6.9) 



in which 



U * " 2 2 K 



(6.10) 



The resulting U2, see (5.7), contains now 9 triplets with circular symmetry, altogether 27 



terms linear in the curvature tensor: 



kU? = 3 



3 (a^R^*? 1 + R p ^A u p J^ - RP^A^J^ 



(6.11) 



The parenthesis represent circular permutations as in (4.6), which explains all those 3's 



If we separate G° into Ricci and scalar curvature parts there are 11 instead of 9 triplets 
and altogether 33 terms. An application of this complicated formula is given in the next 
subsection. 

(Hi) Einstein-Gauss-Bonnet theories of gravity with maximally symmetric boundaries 

In these theories of gravitation, the Lagrangian density is a linear combination of the 
Einstein-Hilbert Lagrangian C\ (a Lovelock Lagrangian of order 1) and the Gauss-Bonnet 
Lagrangian £2 (a Lovelock Lagrangian of order 2) 



£2' = a\C\ + 02^2, 



(6.12) 



in which ct\ and a 2 are self-coupling constants to be found experimentally. The corresponding 
superpotential is thus 

Up = ai E/f + a 2 U 2 av . (6.13) 

In past applications backgrounds have mostly been flat or anti-de Sitter, i.e. maximally 
symmetric spacetimes. In such spacetimes : 



R 



2R 



vpcr 



N(N 
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and R 



R. 

pv - -^g 



1 \IV 



(6.14) 



This makes U2 considerably simpler: 

K- l ";w-* Rt>r. (6, 5 ) 

We have applied this result to calculate the mass of the Einstein-Gauss-Bonnet black holes 
found by Deser and Tekin \7\ with a background curvature 

1 ,_ \ , 2 (N- 1VN-2) 

R/iupa = p {g^pQya ~ 9 ^9 up) where * = ^ > A<0 - ( 6 - 16 ) 

in the case considered in [7], &x = 1 and 

-2l 2 , N 

° 2 = (JV-3)(JV-4)' (6 - 17) 

one obtains the surprising result that 

U^ = -U^\ (6.18) 

But this is indeed the correct result and it gives the positive mass found in Deser and Tekin 
by a quite different but reliable method. 

(iv) Lovelock Gravity in any number of dimensions 

The Lovelock Lagrangian density of order n in N dimensions is as follows: 

. n=p—l 

r — V~9 \ tV\V%—V% v DMl/i2 R M2p-lM2p In 1Q\ 



2 n K 

n=0 



with 



1 < n < p-l where p = \N (for N even) , p = §(iV+l) (for N odd) and 5%£'-%% P = det ( S )• 

(6.20) 

In 1985 Zwiebach [37J speculated that the Lagrangian that would ultimately describe gravita- 
tional interactions and more fundamentally string interactions should have an infinite number 
of terms. Each U„ would give a finite contribution. In particular for spacetimes that are 
asymptotically maximally symmetric, 

f T ,u _ „ (jV-2n)[jV-(2n-l)]---(iV-4)(jV-3) lf 

U n ~ n [N(N-l)]"-i K Ul ' [b ' 

and for ./V — > 00, the factor containing ./V is equal to one. Thus, if the Lagrangian density 

00 

£00 = y^Qn^n, (6.22) 
n=0 
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the superpotential is 

Ug = (jrna n R n ) E7f . (6.23) 



=o 



This sum may converge for suitable as. 

7 Some comments 



a] Regarding U^ a as given by equation (6.5). On a flat background, in Minkowski coordi- 
nates and with £ the Killing vectors of spacetime translations, U^ a , as we observed, is exactly 
the superpotential that Freud [35] found, more than 70 years ago, to calculate mass-energy 
and total linear momentum in Einstein's theory. One wonders why Freud did not calculate 
the angular momentum on the same occasion. 

b] We emphasized several times that the superpotential holds at null as well as at spatial 
infinity and that at null infinity gives the Bondi mass loss in four dimensions [21 J as well as 
in 5 dimensions as we show in Appendix A. We are not aware of a published work on radiating 
fields in Einstein-Gauss-Bonnet theories. Uy is, of course, the superpotential appropriate 
for calculating the Bondi mass loss in that case. 

c] We imposed Dirichlet boundary conditions and found the KBL superpotential in N di- 
mensional Einstein's gravity. Julia and Silva showed that by imposing Neumann boundary 



conditions we are led to Komar's superpotentia 22 Neumann boundary conditions have been 
considered in recent works. See for instance a paper by Kofinas and Olea |24j on Lovelock 
anti-de Sitter gravity. We presume that imposing Neumann boundary conditions would have 
led to Noether charges as considered in Iyer and Wald |13j . 

d] Jacobson and Myers [15] used the first law of thermodynamics to define mass-energy of 
Lovelock black holes. Mass-energy has, however, little to do with thermodynamics because 
asymptotic spacetimes ignore the source of gravity. A direct calculation of the mass-energy 
of Lovelock black-holes, independent of thermodynamic considerations, is naturally provided 
by the superpotential U pa . An approach, similar to that of Jacobson and Myers, was used 
by Gibbons et al [9J to find the mass of Kerr-anti-de Sitter black holes. Needless to say the 
same mass-energy is obtained with the KBL superpotential l^* 7 , see [6]. 

e] We dealt with Lovelock Lagrangians whose variational derivatives have no derivatives of 
the curvature tensor. Conservation laws have, however, been considered for non-Lovelock 



22 A generalization of Komar's superpotential, that does not need a background, to any type of boundary 
conditions, can be found in several papers of Obukhov and Rubilar's ; see for instance [28]. The role of 
backgrounds is here replaced by another ingredient, "generalized" Lie derivatives. 
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Lagrangians by Deser and Tenkin [7\ at spatial infinity and by Wald and Zoupas [36] at null 
infinity. A covariant generalization of Regge and Teitelboim principle to a Palatini or affine 
gravity formulation needs yet to be worked out but do not appear to be trivial unless one 
introduces auxiliary fields in order to reduce second order equations to first order one which 
is always possible while preserving the gauge symmetries, see |17j . 
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Appendix 

A The Bondi loss of mass in 5 dimensional Einstein gravity 



Tanabe, Tanahishi and Shiromizu [33 J calculated the asymptotic solution of Einstein's equa- 

tions in five dimensions at null infinity. They used coordinates similar to those of Bondi 
and calculated, among other things, the Bondi loss of mass-energy. We have calculated the 



use their notations up to minor changes with their signature of the metric and give only the 
quantities needed for the calculation. The Bondi coordinates are x^ = {u, r, x A } thus with 
x A = {9, <p,ij)}, A,B, - ■ ■ = 2, 3, 4 which are the coordinates of a unit sphere of dimension 3, 
£3. In these coordinates, the line element may be written 



ds 2 = 2- + r 2 h AB U A U B du 2 -2e B dudr + r 2 h AB {dx A + U A du){dx B + U B du). (A.l) 



V, B, U depend on all coordinates. The asymptotic expansion in powers of r is as follows: 



V = r 2 + V l r 1 l 2 -£ + 0{r- 1 / 2 ) , B = + 0(r" 4 ) , U A = ^ + 0(r~ 3 ), (A.2) 




mass loss from our formula (6.4) for Ui and have found, of course, the same result. Here we 
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Vi,Bi, U A and Kab are functions of {u, x A } only: 

1 



sin 29 



d e (sin 29 ufj + d^Uf + fyuf 



Bi = -l {C 2 U + CnC 21 + C? 2 + ift + L> 2 2 + D 2 3 ) , 



6 _ 2 
5 



(C n - C 



12, 



, <9</,L> 



C/f sin^ : 
C/f cos0 



tan# 

Dl2 



- (2C n + C 12 ) tan + d e C n + ^ + 



V^12 

sin 9 cos # 



tan# 
1 

tan# 



tan 6* + Oe-Dn + . - + 



2tan0 +a e L>i2 + 



sin 9 cos t 

d^Dn _ ay,(Cii + Ci 2 ; 
sin# 



COS# 



h A B = h AB + ^ +0(r' 5 / 2 ) with /uu = diag{l,sin 2 #,cos 2 #}, (A.3) 



and 



qee = C 11 , qe (p = sin9D 11 , q Sl p = cos 6 D 12 , 

q w = sin 2 6 C 12 , = sin cos 9 D 13 , q^ = cos 2 9 C 13 ; (A. 4) 

The three C's and three D's are functions of {u,x A }, they are not independent: 

C n + C 2 i + C31 = , C12 + C 22 + C32 = , C13 + C 23 + C33 = D 2 U + D 2 21 + (A.5) 

There are thus 5 independent functions associated with the independent degrees of freedom 
of the gravitational waves in five dimensions. £(u, x A ) is an integration "constant". 
The background at infinity is flat. Its line element in our coordinates is 



ds 2 = g^dx^dx" = -du 2 -2dudr+r 2 (d9 2 + sin 2 9dip 2 + cos 2 9di\) 2 ) 



^g = \ r 3 sin 29. 

(A.6) 



To calculate the mass loss we shall need the inverse components of the metric g^ v \ the 
non-zero ones are as follows: 



— 1 + §+^- 7/2 ) - 9 U = ^-i-^ + 0(r-V 2 ) , g^ = ^ + ^ + 0(r-^ 2 ), 



V^_^_Bi 

y»3/2 y»2 ^>3 



^5/2 y7/2 



» 22 - h - % - St) + h ^ - Ca + D - + D 'M + °< r ~ 1I/2 >> 

9 J, = -^—„ ( % + + T^-fl [-C13-D,, + (D 12 - l)Aa]+ 0(r-»'2), 

Ci 2j Di 2 - Z? 23 + D n D 13 ) + 0(r- n / 2 ), 



r J sin C/ \ r 



J Dl2 + ^22_\ + 1 



,33. 



r 2 sin 2 1 



1-- 



,1/2 1 r 3/2 

Ci 2 c 22 



r 5 sin ( 



»3/2 ^5^2 



) + "T 1 ^ (5 ^12 " ^23 + D u + D 2 3 ) + 0(r- 11 /2 ) 
/ r° sin 9 
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- ^ + ^ « " <*• + + + °^ 11/2 )- < A - 7 ) 

gr'"' appears in the Christoffel symbols and in the end we need A A ^ V = T x ^ u — T X ^ U . 
The background calculations are not particularly difficult. One will not need all the A's for 
the energy associated with the time translation Killing vector £ = {1,0,0,0,0}. Moreover 
\J—~g oc r 3 and most of the yj— gUi's go to zero when r — > oo; in the end the only remaining 
term in Ui that contributes to the energy is 

y / —gA 1 0l = \ sin 20 £ + non contributing parts. (A. 8) 

From this follows that the charge 

E = j> £dS 3 (A.9) 

Finally from Einstein's equations R uu = and R ux a = we obtain, like Tanabe, Tanahashi 
and Shiromizu in |33j the Bondi mass loss; with E = Mc 2 

^ = j> [\ [(fluCn + i d u C 21 ) 2 + \ (d u C 21 ) 2 + (d u D n ) 2 + (d u D 12 ) 2 + (d u D n ) 2 ] }dS 3 . 

(A.10) 
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